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A STUDY OF SECOND-ORDER SUPERSONIC FLOW THEORY 


Bv Milton D. -Van Dtke 1 


SUMMARY 

An attempt is made to develop a second approximation to the 
solution of problems of supersonic flow which can be solved by 
existing first-order theory. The method of attack adopted is an 
Iteration process using the linearized solution as the first step. 

For plane flow it is found that a particular integral of the 
iteration equation can be written down at once in terms of the 
first-order solution. The second-order problem is thereby 
reduced to an equivalent first-order problem and can be readily 
xolved. At the surface of an isolated body, the solution reduces 
to the well-known result of Busemann. The plane case is 
considered in some detail insofar as it gives insight into the 
nature of the iteration process. 

Again, for axially symmetric flow the problem is reduced to a 
first-order problem by the discovery of a particular integral. 
For smooth bodies, the second-order solution can then be cal- 
culated by the method of con Kdrmdn and Moore. Bodies 
with corners are also treated by a slight modification of the 
method. The second-order solution for cones represents a 
considerable improvement over the linearized result. Second- 
order theory also agrees well with several solutions for other 
bodies of revolution calculated by the numerical method of 
characteristics. 

For full three-dimensional flow, only a partial particular 
Integral has been found. As an example of a more general 
problem, the solution is derived for an inclined cone. The 
possibility of treating other inclined bodies of revolution and 
three-dimensional wings is discussed briefly. 

INTRODUCTION 

As the linearized theory of supersonic flow approaches full 
development, the question arises whether more exact approx- 
imations are practical. If viscous effects are large, refine- 
ment of the perfect-fluid solution is impractical. If viscosity 
is negligible, however, higher approximations are known to 
yield a closer approach to reality. In intermediate cases, 
an improved solution is desirable in order to assess the 
relative effects of viscosity and nonlinearity. 

The prototype of a higher-order solution for supersonic 
flow is Busemann’s series for the surface pressure in plane 
flow past an isolated body. This simple result is of con- 
siderable value in analyzing supersonic airfoil sections. 
Two terms of the series prove sufficient for almost all require- 
ments; the extension to third and fourth order is chiefly of 
academic interest. 

The aim of the present study is, therefore, to find a second 
approximation, analogous to Busemann's result, for super- 


sonic flow past bodies which can be treated by existing 
first-order theorj’. The natural method of attack, and 
apparently the only practical one, is by means of an iteration 
process, taking the usual linearized result as the first step. 
Several writers have applied this procedure to two-dimen- 
sional subsonic flow. In supersonic flow, as usual, the 
solution is simpler, so that more general problems can be 
solved. 

This paper is based upon a thesis for the degree of doctor 
of philosophy in aeronautics written at the California 
Institute of Technology under a National Research Council 
predoctoral fellowship and under the guidance of Prof. P. A. 
Lagerstrom (reference 1). It was published in revised form 
as NACA TN 2200, 1951 (reference 2). The present version 
has been further slightly revised, in particular to include ref- 
erences to the recent literature. 

ITERATION PROCEDURE 

BASIC ASSUMPTIONS 

The problem to be considered is that of steady three- 
dimensional supersonic flow of a polytropic gas past one or 
more slender bodies. As indicated in figure 1, the bodies 




are assumed either to be pointed or to extend upstream 
indefinitely as cylinders parallel to the free-stream direction. 
In either case, the origin of coordinates can be chosen so 
that all variations in body shape are confined to the half- 
space x>0. Wind axes are introduced, so that for x<0 the 
flow is uniform and parallel to the x axis, with the free-stream 
velocity U and Mach number M. (For definitions of all 
symbols, see the appendix.) 

The bodies are slender, which means that at any point the 
component of U normal to the surface is small compared 


i Supersedes NACA TX 2200, “A Study of Second-Order Supersonic- Flow Theory,” by Milton D. Van Dyke, IS51. 
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with U itself. The symbol « will be used throughout as a 
measure of this smallness. Thus the ordinates of a body 
will be written as e times a function of order unity. Used 
in this way, e serves to distinguish terms of various orders. 

The aim of this investigation is to find a second approxima- 
tion to the solution of problems for which the first-order 
solution is available. The first-order, or linearized, solution 
is defined as the result of keeping only linear perturbation 
terms in the equation of motion. Similarly, the second-order 
solution is the result of retaining second-degree terms ,in 
perturbation quantities. In addition, however, certain of the 
triple products are in some cases found to be as important 
as one or more double products and are therefore also re- 
tained in the equation. It may be noted that the second- 
order solution will not generally consist simply of terms of 
order e and e 2 , though this is the case for plane flow. For 
example, the second-order solution for flow past a body of 
revolution will be found to contain terms as high as 

The flow is assumed to be irrotational and isentropic. 
This assumption is justified in the first- and second-order 
solutions, since the resulting error is found to be at most of the 
order of terms neglected elsewhere. 

EXACT PERTURBATION EQUATION 

Under the assumption of irrotational flow, there exists a 
velocity potential ft. In Cartesian coordinates, the equation 
of motion is (reference 3, equation (39)). 

(c s -fi/)ft«+ (c 8 -iy)ft„+ (c*-ft, 2 )ft„- 

2ftyftjftyg 2ft f G r ft ZZ , 2ftxftyftxy === 0 (1) 

Here the local speed of sound c is related to c„, its value in 
the uniform stream, by 

(Q*+ty+Q*-lP) ( 2 ) 

where y is the adiabatic exponent. The subscript notation 
is used to indicate differentiation. 

A perturbation potential $ is now introduced in the usual 
way. For convenience, however, $ is normalized through 
division by the free-stream velocity so that 

ft= U(x-\-$) (3) 

The perturbation velocity at any point is then the gradient 
of <f> multiplied by U. 

Introducing the perturbation potential into the equation 
of motion gives, after some manipulation, 

p— (2$*+$ I 2 +<v+$.mx+ 

+ 24>,(1 + $*)$«+ 

2(l+*,)*,*, r ] (4) 

where (3=yilf 2 — 1. 

SOLUTION BY ITERATION 

The exact perturbation equation (equation (4)) is com- 
pletely equivalent to the original nonlinear potential equa- 
tion (equation (1)). Simplifying assumptions must there- 


fore be introduced in order to solve it. If it is assumed 
that squares and products of the derivatives of <f> can be 
neglected, the right-hand side of equation (4) disappears, 
leaving the wave equation 

(5) 

This equation is the basis of the linearized or fust-order 
perturbation theory, and its solution is designated by 4> (1) . 

More exact solution of equation (4) by means of iteration 
was first suggested by Prandtl (reference 4). The procedure 
has been applied to plane subsonic flow by Gortler (refer- 
ence 5), Hantsche and Wendt (references 6 and 7), Imai and 
Oyama (references 8 to 10), and Kaplan (references 11 to 
13). Schmieden and Kawalki (reference 14) applied the pro- 
cedure to subsonic flow past an ellipsoid of revolution. Most 
of these writers have considered the stream function rather 
than the potential, which restricts the method to plane or 
axially symmetric flows. The procedure is clearly described 
by Sauer (reference 3, p. 140) for the case of plane floiv. 

The linearized solution $ (l) , subject to proper boundary 
conditions, is taken as the first approximation. Substi- 
tuting this known solution into the right-hand side of 
equation (4) gives 

^J»=F l {x,y,z) ( 6 ) 

where F\ is a known function of the independent variables. 
This is again a linear equation, the nonhomogeneous wave 
equation. A second-order solution subject to, proper 
boundary conditions, can be sought by standard methods. 
The procedure can be repeated by substituting <i>° !) into the 
right-hand side of equation (4) and solving again. Con- 
tinuing this process yields a sequence of solutions <P (#> 
which, under proper conditions, perhaps converges to the 
exact solution. 

A significant feature of this procedure is that in each step 
the left-hand side of the iteration equation js the same. As 
a consequence, the characteristic curves of each iteration 
equation are just the Mach lines of the undisturbed flow'. 
However, in actuality the local Mach lines are usually 
neither straight nor parallel; that is^ the characteristics of 
the original nonlinear equation are curved, in a manner which 
is initially unknown and which depends upon the shape of 
the body. Because of the fundamental role played by the 
characteristics in the theory of hyperbolic equations (see, 
e. g., reference 15, ch. 5; reference 16, eh. 2), it might be 
anticipated that an iteration procedure should be chosen 
such that in each step the approximate characteristics 
would be successively revised so as to approach the actual 
characteristics. For purely subsonic flow', the counterpart 
of such a procedure is knowm to converge under proper 
conditions (reference 15, p. 288-289). Convergence might 
reasonably' be anticipated also in the case of purely' super- 
sonic flow. 

Unfortunately, an iteration procedure in which the 
approximate characteristics arc successively' revised would, 
except in the first step, involve equations with nonconstant 
coefficients. This would greatly, complicate the procedure. 
Fortunately, it will be found that the scheme adopted here, 
which makes no provision for such revision, nevertheless 
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gives an improved solution nearly everywhere in the flow- 
field. 

SECOND-ORDER ITERATION EQUATION 

Henceforth, only the first two steps of the iteration process 
will be considered in detail. In order to eliminate cumber- 
some superscripts, it is therefore convenient to introduce a 
simplified notation for the first two approximations. The 
first-order perturbation potential will hereafter be denoted 
by <p, and the second-order potential 4 by tj> 




(~) 


Introducing these quantities into the exact perturbation 
equation (equation (4)) gives the following second-order 
iteration equation for <f>: 


4>,,+ 0 4 <£.w=-ll a 1 (2vJi+<Pr I + 'P» 1 + <Pt 1 )(<(>xx+<Pi„+ 

<fizz) +'2 <Pi<Piz-!c <Px i pxi+<P, 1 <Pn J r vffzz + 

+ + <Px) <£>it+2 (1 + (.8) 


Since <p satisfies equation (5), the term in 

the right-hand side of equation ( 8 ) can be replaced by 
and the equation for <f> becomes 

, + 4>ti — P*4>tx = -H 2 (2 tp x + ip* + <p* + <P*) Vxx + 

-'PiVxx + Vx'Pxx + <P,*<P V , + VzVxz + 

2<p,(l +pj)^«+2(l +«5 x ) i P»9'i»J 

(9) 


Here the right-hand side contains not only double products 
of perturbation quantities but also triple products. The 
latter can be omitted for plane flow, since they contribute 
terms of smaller order (equal to those found in the next 
iteration). Otherwise, certain of the triple products should 
be retained because their contribution is as great as that 
of one or more of the double products and greater than any 
eontribution from a third approximation. It will be seen 
later that triple products should be retained if they involve 
only derivatives normal to the free stream. Triple products 
which involve x derivatives can be neglected, so that the 
iteration equation becomes 

+ ^ x x=-U 2 { [2+ (7— l)H i l<Px<Pxx + 2<Pt<Pz,+ 

2?is»ic + <Py i ‘P,z J r^Vy<Pi‘Plt J r <PiV2x } (.10) 

Here the triple products which may be important are the 
last three terms on the right-hand size. 


The adiabatic exponent 7 will be found to occur always in 
the combination 


v_ (7 + 1)M 2 


(ID 


i Here + Is regarded as being the complete second-order potential rather than (as fn ref- 

erences 1 and 2) a small correction to be added to the first -order solntfon. 


Introducing this expression in place of 7 gives the final form 
of the second-order iteration equation: 


<t> n +<t> ix— & i 4>xx= -1P[2 {N— 1) fPtpxtpxx + 2 <p l <pxr+2t(>tcp ix + 

i P* < ?ZX ( 12 ) 

ITERATION EQUATION IN OTHER COORDINATES 

In cylindrical coordinates, equation (9) becomes 

t rr +^+^-^ a =w[2iN-l)p i VxVa +2'p rV „+2 ^ „ rf + 

T T L " 


I 1 n 't'Tfm 
(pr <Ptt-T & —J- <Prl r i 


br<pn+ 


f,( -> 2 <P» 1 <Px<PI Y] 

u( >Px~<Pzx’ <Pr <Pxx> ptVxx’ -ps-pxi’ VxVrVxr ) 


(13) 


The terms whose form is indicated at the end of the equation 
are the triple products which will be found to be negligible. 

For conical flows it is convenient to introduce nonorthogo- 
nal conical coordinates (x, t, 6) where 



(14) 


If the body itself is conical, the perturbation potentials are 
reduced to functions of two variables by introducing conical 
perturbation potentials (reference 17) so that 

pCr, t, B )~ *?(£, 0) (15) 

with a corresponding definition for <j5. The derivatives are 
given by 


¥>*=?— fat 

f 4 _ 

¥>«=- <Pu 

* 

9=r« = j3p«l 




P- 

<Prr — x <Pu 

4>i# 

►> 

(16) 

9=1 — X ft 

<P»=x<Pu 



Si* 

1 

II 

sS 




with the same relations connecting 4> and <j>. The iteration 
equation (equation (13)) becomes 


(1 -f 2 ) 5,,+^+^= IP [ 2 (A r — 1 )P(v-tj c ) 2 


{ 

j- 5i(pi — t<p,i) J rp r <pi' !: pit-\-2 <p t <p»<pt > — 

pt Vi<p»*+p. f <pi r <pt»+0^J i (<p— fpt)V«o 


PvFptt, ¥»#Vti! ji (<P~ 




(17) 


Here the grouping of terms corresponds to that in equa- 
tion (13). 
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BOUNDARY CONDITIONS 

Physical considerations suggest that the flow should satisfy 
the following conditions: 

Tangency condition: The resultant flow is tangent to the 
surface of the body. 

Upstream condition: All flow perturbations vanish every- 
where upstream of the body. 

Those two requirements are sufficient to determine the 
solution. The first imposes one condition along the time- 
like surface of the body, and the second may be regarded as 
imposing two conditions on a spacelike surface. This is the 
case of mixed boundary conditions (reference 15, p. 172) 
and leads to a determinate solution (see reference 16, p. 85). 

The tangency condition may be written formally as 

grad Q - grad S=0 at £=0 (18) 

where S= 0 is the equation of the surface of the body. In 
a more useful form it becomes, for the first- and second-order 
problems, 

ip e = (slope) (l-f-pi) (on the surface) (19a) 

<j> c — (slope) (l+ff’r) (on the surface) (19b) 

Here <#> £ means the cross-wind component of the normal 
derivative of <t> at the surface of the body. In plane flow 

is 4 >„, and in axially symmetric flow is <f» r . The slope of 
the body is measured with respect to the free-stream 
direction. 

In first-order theory, the tangency condition is frequently 
approximated by neglecting <p x compared with 1 in equation 
(19a), which causes only a second-order error. Correspond- 
ingly, in the second-order tangency condition (equation 
(19b)), <t>z can be replaced by its first-order counterpart <p x 
with only a third-order error. Thus the tangency condition 
simplifies to 


The term containing <j>, in equation (21b) accounts for the 
fact that to second order the cross-wind component <t> e may 
not be vertical. Corresponding simplifications can be made 
for any quasi-cylindrical body, which is a body whose entire 
surface lies near to a cylinder parallel to the free stream. 

Fmalty, the approximation which led to equations (20) 
may be adopted in addition to that just discussed for planar 
systems, in which case the tangency condition simplifies 
further to 

<p y — Y x at y = 0 (22a) 

4>y= Ti(l + ^x) — Y<pyy-\- at y=0 (22b) 

The upstream condition implies that both <I> and <P X vanish 
at the plane z=0, which completes the boundary conditions 
required. For the first- and second-order problems the 
upstream condition is therefore 

<p=tp z = 0 at /=0 (23a) 

<#>=<£*=0 at r=0 (23b) 

DETERMINATION OF PRESSURE 

When the potential field has been determined, the net 
velocity q at any point is given by 

g* x= (ZJ+uy+v l +w* (24) 

where 



y > c = slope (on the surface) (20a) 

<j> e = (slope) (l + «p z ) (on the surface) (20b) 

This approximation will not be made except for plane flow. 

Another approximation in the tangency condition can be 
made for planar bodies. A planar body is one whose entire 
surface lies near to a plane, parallel to the free stream, say 
the plane y=0 (reference 17, p. 52). Thin flat wings are 
planar bodies, whereas slender pointed bodies of revolution 
are not. For a planar body the first-order tangency condi- 
tion can be imposed at the plane rather than on the surface 
of the body. Correspondingly, the second-order tangency 
condition can also be imposed at the plane provided tliat the 
difference in <#>„ between the plane and the surface of the 
body is accounted for approximately by retaining the second 
term in its Taylor series expansion about y=0. Further- 
more, to first order <p t is <p„ for a planar body. Tims if the 
surface of a planar body is given by y= Y(x,z), the simplified 
tangency condition is 

<p a =Y x ( l + <p x ) at y=0 (21a) 

4> y -hY<f>yy— F,<#> t =Fx(l+<fce) at y=0 (21b) 


in Cartesian and cylindrical coordinates, respectively. Be- 
cause the flow is assumed to be isentropic, the pressure 
coefficient is given by 

c * m r‘ P “ = ^{[ 1+L 2 _MJ ( 1 ”F i )] -1 S (20) 

2 P, ° U 


where p„ and are the free-stream pressure and density. 

It is the practice in linearized theory also to simplify the 
pressure relation. Substituting equation (24) into equation 
(26) and expanding gives 



t>*+w ! , at u 2 , p s +w ! , M Y v*H-w*Y , 

U% TP u jji 4 \ U* J + 


v? u 2 v 2 +w 2 u /ri+wV /'v i +w 2 y~ 

u Lu s ’ u* u 1 ’ u\ u 2 )’\ u 2 ) _ 


(27) 


All the terms shown here explicitly may give contributions of 
second order; the remaining terms whose, form is merely 
indicated are always of higher order. In linearized theory 
only the first term is ordinarily retained. This is satisfactory 
for plane flow or flow past planar systems, since the contribu- 
tion of the remaining terms is definitely of higher order. 
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In fact, for plane flow past an isolated body it happens that 
the next two terms cancel identically. 

For slender bodies such as a cone, however, orders of 
magnitude are not so clearly distinguished. Busemann 
suggests (reference 18) that the second term, (p'-hu^/t 71 , 
is then sufficiently large compared with the first that it 
should be retained also. This view is supported by Lighthill 
(reference 19), who shows that the resulting solution is 
correct up to the order of the quantities contributed by the 
second term. Again, the third term, phi*/!.*, which is also 
the square of a perturbation quantity, is comparable with 
the second at high Mach numbers (where LightliilTs order 
estimates become invalid) and might logically be retained. 
Having gone this far, it may be simpler to use the exact 
relation. 

Each of these four possibilities for the first-order flow past 
a 5° cone is compared with the exact solution (reference 20) 
in figure 2. The series (equation (27)) is seen to alternate in 
this case. It converges so slowly, however, that at moderate 
Mach numbers (say, near d/=\2) where first-order theory 
is most accurate, linearizing the pressure relation introduces 
much greater errors than linearizing only the equation of 
motion. Adding each of the quadratic terms in turn causes 
changes nearly as great as the error due directly to non- 
linearity in the equation. 

The point of view to be adopted here is that calculating 
the velocities and calculating the pressure are two essentially 


distinct operations. A certain degree of approximation may 
be necessary in order to solve for the velocities, but the 
pressure relation need not then be approximated to the same 
extent simply for the sake of consistency. For it may happen 
that the resulting errors (though of the same mathematical 
order) are greater than those due to the original approxima- 
tion. Indeed, this is evidently the case at moderate Mach 
numbers in the first-order solution for a cone and will be 
found true to a greater extent in the second-order solution. 
Moreover, in the second-order solution so many terms of 
equation (27) must be retained that it is usually simpler to 
use the exact relation. For these reasons, the exact pressure 
equation (equation (26)) will be used throughout except in 
the case of plane flow. 

EOLE OF A PABTICULAB LNTECBAL 

The second-order iteration equation can be attacked by 
standard methods, and in the case of plane flow a solution 
can be found directly. For plane and axially symmetric 
flows, however, it will be found that a particular integral of 
the iteration equation can be written down at once in terms 
of the first-order solution. This solves the problem because 
the complete solution consists of a particular integral plus a 
solution of the homogeneous equation, and the latter can be 
obtained by existing methods. That is, the second-order 
potential may be written as 

<t>=f+X (28) 

where 
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\p is any particular integral of the nonhomogeneous iteration 
equation 

X is the complementary function which serves to correct the 
tangency condition 

The problem for x is the usual first-order problem for which 
methods of solution are assumed to be known. 

The role of the particular integral is to transfer the non- 
homogeneity in the problem from the equation, where it is 
troublesome, to the boundary conditions, where it can be 
handled by existing theory. For linear partial differential 
equations it is always possible in principle to transfer non- 
homogeneities in this way from the equation to the boundary 
conditions, and vice versa, by adding a suitable function to 
the dependent variable (see reference 21, p. 236). 

Since the particular solution 4* will be found in terms of the 
first-order solution, it will vanish upstream of the plane a?=0. 
Then the complementary function must also vanish there, so 
that the upstream condition for x is 

x(0,3/,z) = x*(0,Jt,2)=0 (29) 

From equation (19b), the tangency condition for x is found 
to be 

tc+Xc— (slope) (l+£r+Xc) (on the surface) (30) 

or, in the case of a planar body given by y= Y(x, z), from equa- 
tion (22b), 

4'i/ J rXt= F r (l + ^x) — Y<p„+Y g <p, e.ty=0 (31) 

1 1 should be noted that, although <j> is of the same magni- 
tude as <p, this is not necessarily true of either ^ or x alone. 

PLANE FLOW 

The second-order solution for conditions at an isolated sur- 
face in plane supersonic flow was given by Busemann (refer- 
ences 22 and 23). By using the iteration procedure, the 
solution will now be found throughout the flow field, includ- 
ing the case when several bodies interact. 

The solution for plane flow is of interest chiefly insofar as 
it serves as a guide in more complicated problems. In par- 
ticular, it provides insight into such details of the iteration 
process as the question of its success and the effect of sharp 
corners. 

PARTICULAR INTEGRAL FOR PLANE FLOW 

The first-order equation for plane flow is 

— 0V*r=O (32) 


To tliis must be added a solution x of the homogeneous 
equation (equation (32)), which has the form 

X=H(x-py)+J(x+fly) (36) 

where II and J are functions determined by the second-order 
boundary condition. 

For flow past a single boundary (such as one surface of an 
airfoil) the first-order potential (equation (33)) contains 
only one or the other of the functions h and j. In this case 
<Pv<Piv= P 2 <Px<Pzx so that the iteration equation reduces to 

<j>n—p 2 <t>xx=2M 2 f3 ! Ntpx<f> tx (37) 

The particular integral may then be simplified to 

4'= M 2 j y<Px<p, (3 8) 


and the complementary function contains only II or J, 
according as the first-order solution contains only h or j. 

FLOW PAST A CURVED WALL 

As an example of the application of the particular solution, 
consider flow past a wall which at some point begins to 
deviate slightly from a plane (see fig. 3). The wall can be 
represented by 

y=Y(x)=eg(z) (39) 


where e is a parameter small compared with unity and g(x) 
is a continuous function of order unity which vanishes for 
zgO. 

This is a planar body, so that the tangency condition is 
given by equation (22a). Consequently, the first-order 
problem is 


l Pw—P I <Pxx=0 ' 

<P v (x,0)=tg\x) 
<pt0,y)=<Px(0,y) = 0. 


(40) 


where the prime indicates differentiation of g with respect to 
its argument. The solution is 


<p=—j;g( x -Py) (41) 

Substituting this first-order solution into the right-hand 
side of equation (37) gives the iteration equation 


The general solution is 

<p(x,y)=h(x— Py)+j(x+$y) (33) 

where h and j are functions determined by the first-order 
boundary conditions. 

In the iteration equation, all triple products can be neg- 
lected, and equation (12) becomes 

[(AT— 1) j3 : ViP«+ Wzy] (34) 

It can be verified directly that a particular integral of this 
equation is given in terms of the first-order solution by 

4>= M 2 <pz 1 ~y) v y<p *\ ( 36 ) 


P 2 4>xx- 2M 2 Ne‘g'(x — Py)g' ' (x - py) (42) 


y 
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According to equations (38) and (36), the solution is 

l/ 2 V 

<p=$+x= — ^y[g\x— /3y)l*+H(jr— j8y) (43) 

Imposing the approximate second-order planar tangency 
condition (equation (22b)) gives 

[* / (*)l f -«(*)« /, (i)| (44) 

so that 

//(jr)=-«‘|«(x)*'(x)+ J/ ^pS i £ [g'i&m] (45) 

The complete second-order perturbation potential is therefore 
)-e 2 |f(x-/Sy)«'(x— /Jy)+ 

ytf X*-m a + M&IW } ( 46 > 

The same result can be found by solving equation (42) 
directly, using the impulse method (reference 15, p. 164). 

Flow quantities at the surface of the wall can be related 
to their values at the plane y = 0 by expanding in Taylor 
series and discarding terms in e 3 . In this way the streamwise 
velocity perturbation at the wall is found to be given by 

„ < I f2\r_9 

- - 2p* - r ' ( 47 ) 

The pressure coefficient at the wall can now be calculated 
from equation (27) which, upon replacing A r by its value from 
equation (11), gives 




(48) 


This means that to first order the actual characteristics are 
approximated by the Mach lines of the undisturbed flow. 

For the second-order solution, a closer approximation to 
the characteristics could easilj - be found. It can be shown 
that if the first-order streamwise perturbation velocity at any 
point in the flow is ti a> , then the revised local values of 
Mach number and /3 are approximately 3 

M (l) =M [l +0*(xV- 1) (50a) 

|8 (U Af®*- l=/3[l + M*(X- 1) (50b) 

Combining this result with the first-order solution (equation 
(41)) shows that the revised Mach lines have the slope 

(51) 

However, because the iteration scheme adopted does not 
allow for such revision, these are not actually the character- 
istics of the second-order equation. Instead, the character- 
istics of equation (42) are still the original Mach lines of the 
undisturbed flow. 

Physically, the characteristics are lines along which 
discontinuities in velocity derivatives are propagated, and 
this definition is completely equivalent to the mathematical 
one (reference 15, p. 297). Therefore, in the second-order 
solution given above, discontinuities in acceleration will 
occur along the original characteristics rather than, as they 
more properly should, along the revised characteristics. 

Suppose, however, that no such discontinuities occur. 
For flow past an isolated body the downstream character- 
istics are also lines along which the velocity is constant, 
provided that shock waves do not appear. Setting 


This is the well-known result of Busemann (references 22 
and 23). To second order, the surface pressure coefficient 
depends only upon the local slope. 


BOLE OF CHARACTERISTICS 


It was pointed out previously that because of the under- 
lying significance of the characteristics for solutions of 
hyperbolic equations, it might be expected that the approxi- 
mate characteristics of the iteration equation would have to 
be revised successively at each stage. However, an iteration 
procedure was adopted which involves no such revision. 
It is therefore pertinent to inquire in this simple example 
what roles have been played by the actual and approximate 
characteristics. 

The flow past a single curved wall is given (until shock 
waves form) by a simple wave or distributed Prandtl-Meyer 
expansion. Of the tu f o families of characteristics, those of 
primary importance in a simple wave run downstream away 
from the wall. We therefore confine attention to that family. 

For the first-order equation, equation (32), the conven- 
tional theory (e. g., reference 15, ch. 5; reference 16, ch. 2) 
shows that the characteristics of the downstream family 
are the lines of slope 


dx (S 


(49) 


d4>r=<t>xxdx+<t) xll dy=0) 

d <t>, = $x,d x + y = 0 ) 

it is seen that the velocity is constant if 
dy ^jx_ 4>xy 

di <t>n 


(52) 


(53) 


For the second approximation (equation (46)) the velocity 
is constant along lines of slope 


dy_l 
cfx d 


1 




(54) 


which, according to equation (51), are the revised character- 
istics. Consequently, although the characteristics have not 
been revised in the mathematical sense, the solution behaves 
physically as if they had, so long as discontinuities do not 
occur. The question of discontinuities will be considered in 
the next section. 

The connection between the original and revised character- 
istics can be interpreted physically. The right-hand side of 
the Iteration equation may be regarded as representing the 

• As Implied by the superscript, these are regarded as first-order values, because they are 
determined Irom the first-order solution. From this point of view, M and 0 are zero-order 
quantities. 
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effect of a known distribution of supersonic sources through- 
out the flow field. The influence of this source distribution 
spreads downstream along both families of original character- 
istics. The resulting velocity changes are just such that to 
second order the velocities become constant along the 
revised rather than the original characteristics. 

Finally, it is interesting to note that the second-order 
potential is constant on lines which bisect the original and 
revised characteristics. For, setting 

d<p—<l) X dx+<t>ydy=0 (55) 

<j> is found to be constant along lines of slope 

FLOW PAST A CORNER AND A PARABOLIC BEND 

A simple case in which discontinuities may occur is that of 
flow past a sharp comer. The exact solution is known to 
involve an oblique shock wave with attendant velocity dis- 
discontinuities for compression and a continuous Prandtl- 
Meyer fan for expansion. 

Denoting the tangent of the deflection angle by e, positive 
for compression (see fig. 4), the function g(x) appearing in 
equation (39) is 



FiiJOBE 4.— Flow past a comer. 


From equation (46) the second-order perturbation potential 
is found to be 


corresponds to a shock wave, for it is known that an oblique 
shock bisects the Mach lines to a first approximation (refer- 
ence 16, p. 354). In the case of expansion, this typo of 
correction cannot be justified since it would involve con- 
tinuation of the free-stream potential across a true charac- 
teristic. Instead, a Prandtl-Meyer fan must be inserted. 



Fioube 5 .— Mach lines before and after adjustment of potontlal discontinuity. 


Evidently the iteration process is successful except within 
an angular region of order « lying near the Mach line from the 
corner. In particular, the pressure is given correctly every- 
where on the surface of the wall. 

It is enlightening to observe that the alternative method of 
iteration, in which the cliaracteristics are successively re- 
vised, fails in the same region. The potential is double- 
valued over a fan-shaped region in the case of compression 
and is left undefined over a similar region in the case of ex- 
pansion (see fig. 6). The same artificial corrections are 
necessary to complete the solution. 



Figure 6.— Second-order flow past a comer using revised characteristics. 


- e t M 2 N 

0y)+^(z-02 /) — ( * X ( 68 ) 

to the right of the line x=Py and zero to the left. Conse- 
quently, in either compression («>0) or expansion (e<0) 
the second-order potential suffers a discontinuous drop along 
the Mach line from the corner, of strength proportional to 
the distance from the corner. Such a discontinuity cannot 
be admitted, which indicates that the iteration process fails 
in this region. 

In the case of compression, the solution can bo corrected 
by analytically continuing the perturbation potential up- 
stream until it can be joined continuously to the free-stream 
pot ential. (This is permissible since the line of discontinuity 
is not actually a characteristic.) From the result of equation 
(56) the juncture is seen to occur along the line from the 
corner which bisects the upstream and downstream Mach 
directions, as indicated in figure 5. The adjusted juncture 


Consider next flow past a parabolic bend which is repre- 
sented by 

y=~tx 2 x>0 (59) 

From equation (46) the second-order perturbation potential 
is found to be 


>(*— /fy)* 


M 2 (N+l)—3 


'(x-pyf 


e*y{x- py)* 
(GO) 


The potential and also the velocities are continuous, so tliat 
the previous difficulties do not occur. The acceleration is 
discontinuous across the original characteristic x=f)y, which 
in this case happens to be also a revised characteristic. 
However, a new complication arises. It is known that, in 
the exact solution for the compressive case, the character- 
istics form an envelope, as shown in figure 7. Inside the 
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cusp the potential is triple-valued (reference 16, p. Ill), so 
that a shock wave must be inserted. This envelope must 
also arise in the second approximation since the character- 
istics are no longer parallel. However, the second-order 
potential given by equation (60) is single-valued, so that it 
cannot predict the formation of an envelope. Again the 
iteration process fails in a part of the flow field. 

It can be seen that the alternative iteration process, using 
revised characteristics, will produce an envelope. 

CONVERGENCE FOR PLANE FLOW 

The examples just considered indicate that the success of 
the iteration procedure should be carefully investigated. A 
step of an iteration process may be considered successful if, 
in some sense, it significantly improves the solution. In 



particular, one is interested in the success of the second- 
order solution. 

It should be noted that a divergent process may be suc- 
cessful for many steps and that, on the other hand, con- 
vergence does not necessarily imply success. In practice, 
however, one would expect a convergent process to be suc- 
cessful. As used here, success is a subjective notion, not 
amenable to analysis. Consequently, only the convergence 
of the iteration procedure can be considered in any detail. 

Unfortunately, proofs of sufficient conditions for converg- 
ence have not been obtained, even in the case of plane flow. 
However, the preceding examples suggest certain conjectures 
regarding convergence. These will be stated and some argu- 
ments for their plausibility will be advanced. 

For flow past a slightly curved plane wall represented by 
y—eg{x), the solution obtained by iteration using the revised 
characteristics is conjectured to converge in any bounded 
region adjacent to the wall, provided that 
(a) e is sufficiently small, 
lb) g(x) is continuously differentiable. 

If tj(x ) has only a piecewise continuous derivative, the 
convergence holds except possibly in fan-shaped regions 
springing from each corner, which lie near the original 
Mach line and subtend an angle of order e. 

For the iteration process actually adopted, in which the 
characteristics are not revised, the first n steps are conjec- 
tured to form part of a convergent process, provided that 
(a) e is sufficiently small, 

(b') g(x) has continuous derivatives up to (n — 1) st order 
if the potential is required, nth order if the velocities are 
requited. 


If condition (b') is satisfied only piecewise, the result holds 
except possibly in fan-shaped regions springing from each 
comer. 

In the first case, condition (a) is necessary in order to 
insure that the solution be unique, as is clear from the exam- 
ple of the parabolic wall. The preceding examples also 
show that condition (b) is necessary. 

If the sufficiency of these two conditions is assumed, their 
connection with condition (b') in the second case can be 
illustrated by analogy with a mathematical model which 
retains the essential difference between the two iteration 
processes; namely, that the correct characteristics are not 
used in the method actually adopted. 4 Consider the first- 
order problem given by equation (40): 


<Pvt~ c Pzx—0 

< P ,(x,Q)=eg'(x) 

<piO,y)=<p c (0,y)=0 


(61) 


where we have taken /9=1 for convenience. The solution 
(equation (41)) is 

<p=-tg(x—y) (62) 

Now we attempt to solve this problem using characteristics 
which differ from the true characteristics by 0(e). Thus 
we consider the equivalent problem 


<Pn— (1 — «) 

<P,(x,0)=eg'(x) 
<p(0,2/)=s>*(0,i/) = 0 


(63) 


and solve by iteration. In the first approximation the 
right-hand side is neglected, so that the differential equation 
becomes 

— (1 — e)<p x J 1) = Q (64) 

which has the solution, subject to the boundary conditions, 

<P W = — eg{x — Vl — ey) (65) 

Substituting this into the right-hand side of equation (63) 
gives the iteration equation for the second approximation: 


^'*-(l- S W*=- e Y'(x-Vl (66) 


Using the impulse method (reference 15, p. 164) gives the 
solution, subject to the Ixmndary conditions, as 

<p® = — egfx— Vi— ey)-r ^ <hjg' (x—y/l — ty) (67) 


But this is just the Taylor series expansion, correct to 0(£*) r 
of the true solution (equation (62)). Subsequent iterations 
add additional terms to the Taylor series expansion. Hence, 
despite the use of slightly incorrect characteristics, the itera- 
tion process converges to the correct solution. The connec- 
tion between conditions (b) and (b') is thus seen to be that 
the existence of sufficiently many continuous derivatives 
compensates for the fact that the wrong characteristics 
are used. 

* This model was suggested by Prof. C. E. DePrima of the California Institute of 
Technology. 
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AXIALLY SYMMETRIC FLOW 

Before discussing the general solution for bodies of revolu- 
tion, it is convenient to consider tlie simple problem of a 
cone. In this case the second-order solution can be found 
directly. The results will be useful in indicating which 
triple products should be retained in the general case. 

FLOW PAST A CONE 

Consider flow past a slender cone of semivertox angle 
tan -1 e, as shown in figure 8. The flow is conical and axially 
symmetric, so that the iteration equation is given by equation 


and from equations (19b) and (23b) the corresponding 
boundary conditions are 

i 30 / =e(l +4>-t4> > ) at f=/3e (71b) 

4 >=cj >'= 0 at t= oo (71c) 

Equation (71a) can again be solved using the integrating 
factor t/ijl—t*. The various integrals encountered can be 
treated by integrating by parts one or more times. Imposing 
the upstream condition (equation (71c)), the complete 
conical second-order perturbation potential is found to be 



0=— - B (sech -1 1— — Jlsech -1 tf— 

(N- f l)sfT=? sech ~H pA 1+0 [ £ ° (sech- 1 /)»] 

" (72) 

From equation (1G), the stream wise and radial velocity 
perturbations are 


u 


— = — 5 sech -1 f+A. 2 M 2 | (seeb+f)*— 


1 ( 




(73a) 


(17) with 8 derivatives omitted. Including the boundary 
conditions from equations (19a) and (23a), the first-order 
problem is 

(1 -<*>?" (/)+^=0 (68a) 

8v'=f(l +v-tv') at t—fie (68b) 

i?=5>'=0 at t= oo (68c) 

The differential e quatio n can be solved by using the inte- 
grating factor t/Vl— t*. Imposing the upstream condition 
(equation (68c)) gives the well-known result 

£=—.4 (sech- 1 f-Vl-* 2 ) (69) 

which is understood to vanish except within the downstream 
Mach cone, where fgl. The tangency condition (equation 
(68b)) is satisfied by putting 

^ o 

A = Vf—+c 5 + e 2 sech- 1 (ySe) (70) 

At the Mach cone (£=1), all velocity perturbations vanish, 
so that the first-order solution predicts no deflection of the 
shock wave from the Mach cone (see reference 16, p. 403). 

Substituting the first-order solution into the iteration 
(equation (17)) gives 


(l-f)0"(f)+^=;i 2 M J 


[2 (AT- 1) 5^+2 1- 


Vi-f 2 


p*A 


Vi-f* 


] 


(71a) 


A tj Vi — f 2 
8 U~ a t 


+A*AP 






„ Vl — sech -1 1 

~ 2 1 h 

.^141 (?±aiiEE 

/f3f> 4 t 


(73b) 


T^he constant B must be adjusted so as to satisfy the tangency 
condition (equation (71b)). In actual compulation it is 
easier to. adjust B numerically in exactly this fashion rather 
than to calculate it from the cumbersome expression which 
could be written down. The pressure coefficient at any 
point can then be calculated from equation (26). 

The last term in the bracket in equation (71a) is the triple 
product /fl?,*?,, which is retained in the iteration equation 
(equation (17)). Its retention is now’ justified by noting 
that its contribution — the last term in equation (72) — is of 
the same order as other terms near the surface of the cone 
(t=p s). Actually, it contributes a second term, which has 
been neglected since it is at most of order ^sech ~ l @t. It 
can also be verified that the other triple products, whose 
form is indicated at the end of equation (17), are in fact 
negligible since they contribute at most terms of order 
^(sech-'/Je) 2 . Consideration of a further iteration indicates 
that a third approximation would add terms of order 
€ 5 (sech- 1 j3e) 3 , which is greater than the terms just neglected. 

The second-order result for surface pressure coefficient is 
compared in figure 9 with the exact solution (reference 20) 
for cones of 5°, 10°, 15°, and 20° semivertox angles. Also 
shown for comparison are the first-order results based upon 
the exact isentropic expression (equation (26)) for the pres- 




Fiii c: re ti.— Comparison of virions solutions for pressure on a cone 


sure coefficient. The second-order solution is seen to provide 
a much better approximation throughout the range of 
Mach numbers up to the point at which the Mach angle 
equals the cone angle, beyond which the perturbation solu- 
tions have no physical meaning. 

SHOCK-WAVE POSITION FOE CONE 

The solution for plane flow fails near the Mach wave from 
a comer, which suggests that the second-order solution for 
the cone may likewise fail near the Mach cone from the 
vertex. In the plane case, nevertheless, a first approxima- 


tion to the shock-wave position (and hence to the entropy 
change) can be calculated from the velocity perturbations 
near the Mach wave. We now consider whether this is 
true also for the cone. 

It was noted before that first-order theory predicts no 
disturbance at the Mach cone and hence no shock wave. 
According to second-order theory the velocity perturbations 
just behind the Mach cone are (equations (73)) 
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so that the perturbation is normal to the Mach cone. Here 
A (equation (70)) has been approximated by 6 s . From 
equation (50b) the cotangent of the revised Mach angle 
just behind the Mach cone is found to be 

j8 w> =0[1 — 2MW(Af— 1)«*] (75) 

The upward stream inclination there is approximately 
(vlU)t-i, so that the Mach lines have the slope 

3~(l+2 MW) . (76) 


Now if this can be taken to be the slope also of the revised 
characteristics just behind the shock wave, then the slope 
of the shock wave differs from that of the original Mach 
cone by 

1 t M*N* 4 (t+1) 2 M 8 4 

0 6 40 s € 


tan X— 

p 


(77) 


This problem has been treated rigorously by Liglitkill (refer- 
ences 24 and 25) and by Broderick (reference 26), who find 
that actually 

1 .3_ fr+p £,, (78) 


tan X— 


P 8 fl 5 


which is one and one-half times the preceding result. The 
discrepancy means that the second-order solution fails near 
the Mach cone. The nature of this failure and the proper 
method of remedying it liave recently been studied by Light- 
hill (reference 25). 

The entropy increase through a weak oblique shock wave 
is proportional to the cube of its inclination away from the 
Mach lines. Consequently, the entropy rise through the 
shock wave from a cone is 0(« 1S ), as noted by Light-hill 
(reference 24). 


PARTICULAR INTEGRAL FOR AXIALLY SYMMETRIC FLOW 

Consider flow past a body of revolution which is either 
a slender pointed body with nose at the origin or one which 
extends indefinitely upstream with constant radius a for 
sc ^0 (see fig. 10). The latter shape corresponds to the 
external flow past a sharp-edged, open-nosed body with 
supersonic flow at the lip. With slight modification the 
subsequent development can be applied to internal flow as 


r 



well. 

The meridian curve can be represented in 

tlie first 

case by 

r=R(x)=€g(x) x^Q 

(79a) 

and in 

the second by 



r./ \( a x ^° 

T X (a + eg(x) x£0 

(79b) 

Here e is again a parameter small compared with unity, and 
g(z) is a continuous function of order unity. 

The first-order problem is 


<Prr+~— 0V«=O 

(80a) 


V[ (x,R)=R f (x)[l+<p t (x,m 

(80b) 


V>(fi,r)=<pA0,r)=0 

(80c) 

The solution is known to be (reference 27) 



r* nm - 

Jo V(x — g) 2 — / 3*r s 

. ----- 


. . z—b 
/VjOSIl"* - 3 — 

— I 9r /(x-jSr coshu)dit 

(81) 


The second form is useful for carrying out differentiation, 
after which the first form can be restored. The derivatives 
which will be required are 


fcosh-' — - 

<p t —— J 9 j'{x— fir cosh v)du= 


Ri-Bt 

J" 


nm 

V(*-0*— 

(82a) 


•’■"'’I 

if 

rjb 


1—6 

cosh-i 


j r (jc — Pr cosh u ) cosh u du= 


1 r x ~* r {x-Qf(Qdj 
V(z-£)’-0 s r* 


(82b) 


In carrying out the differentiation the fact has been used 
that /(ft) =0 for a body with finite slopo. With coordinates 
as shown in figure 10, the lower limit of integration b is 0 
for the pointed body and —flu for the semi-infinite body. 
The function /(x) may be interpreted physically (aside from 
a numerical factor) as the strength of a supersonic line source 
along the x axis. It is determined by the tangency condition 
(equation (80b)) which gives the following integral equation 
of the Volterra type for/': 


p-* aw (x-e/wt 
Jo Vk-S’-W*) 


R(x)R 




x-BR {*) 


urn. 




] 


(83) 


The second-order iteration equation is found from equa- 
tion (13) to be 

— P 2 (p a = M 2 [2(N— l)/JV*#tt + 2^xr + P r Vrr + 

T 

0(¥J*V*r, <Pz<Pt<Pzt)} 


Figpbi 10.— Flow past bodies of revolution. 


(84) 



A STUDY OF SECOND-ORDER SUPEESONIC FLOW THEORY 


501 


The solution for the cone suggests that the terms indicated, 
at the end of the equation are negligible. 

It will now be shown that a particular integral of this 
equation is given in terms of the first-order solution by 


known as the slender-body approximation (reference 27). 
For slender bodies, the source strength f(x) appearing in 
equation (81) is found to be approximately proportional to 
the rate of change of cross-sectional area. Thus 


^(r,r)=ilp[^+A>^)-^^ 3 ] (85) 

The first group of terms contributes the first two terms on the 
right-hand side of equation (84) , as can be verified by direct 
substitution. The last term in equation (85) accounts for 
the term tp 2 <p„ as follows: 


<Pr 1 <Prr — ^ i <Pr(^>'<p a ip rr + <p r ipa—^r<p z /j 


(86) 


-where repeated use has been made of the fact that <p satisfies 
equation (80a). The last group of terms consists of triple 
products involving x derivatives, which have already been 
neglected in equation (84), so that the result is proved. 

The complementary function x is a solution of equation 
(80) and can be written as 

(87) 


Using equations (82) the second-order tangency condition 
(equation (19b)) is found to be 






V(* 

F'(Qd£ 


Vfcr-tf-jMP. 


0 


( 88 ) 


which is again a Yolterra integral equation. 

METHODS OF SOLVING INTEGRAL EQUATION 

Discovery of a particular integral for bodies of revolution 
reduces the second-order problem to the same form as the 
first-order problem; namely, the solution of a Yolterra in- 
tegral equation. Y arious methods of attacking this problem 
are listed by Hayes (reference 17, p. 140). 

An indirect method consists in assuming that the unknown 
source strengths in equations (81) and (87) can be represented 
by a few terms of a polynomial, for example, 

f(x)=C l x+C^+ . . . +C a x* (89) 


The resulting solutions were introduced in a more formal 
manner by Hayes (reference 17, p. 38), who has discussed 
their properties in detail. The first term alone gives the po- 
tential for the cone, (equation (69)). Additional terms give 
the solution for simple families of shapes. However, the 
method is not suitable for bodies having discontinuities in 
slope or curvature. Consequently, a more direct procedure 
is desirable. 

von K&rmdn first introduced an asymptotic solution of the 
integral equation (equation (83)) which has come to be 



Figure II. — Equivalence of polygonal source and sum of conical sources. 


Lighthill has shown (reference 19) that if R(x) and its first 
two derivatives are of order e and R' is continuous, this de- 
termination of /(x) is correct to the order of terms retained 
in the first-order solution. For purposes of the second-order 
solution, it can be shown that f(x) may be determined in this 
way only if the first four derivatives of R are of order e and 
R" is continuous. This means that the body must have 
continuous curvature, which is a severe limitation. More- 
over, the slender-body approximation is found generally to 
cause unnecessary loss of accuracy even though the mathe- 
matical order estimate of the error is small. Consequently 
this approximation should be avoided if possible. 

The most satisfactory way of solving the integral equations 
is to use a step-by-step numerical procedure. In first-order 
theory the usual method, introduced by von K&rm&n and 
Moore (reference 28) , is to assume that the unknown source 
distribution can he approximated by a polygonal graph. 
This is equivalent to superimposing a number of conical 
source lines of different strengths, each shifted downstream 
with respect to its predecessor, as indicated in figure 11. 
The latter viewpoint is more convenient for computation. 
The strengths of the source lines are determined in succession 
by satisfying the tangency condition at a series of points on 
the surface of the body. The details of this prodedure are 
clearly, explained in reference 3, page 77. 

For purposes of a second-order solution, this procedure 
must be modified in one respect. Unless the source distri- 
bution /(a.) actually has corners, it must not be approximated 
by a polygon. The reason is that a comer corresponds 
locally to adding a conical source line which would, according 
to the solution for the cone, produce false second-order dis- 
continuities in velocity and pressure across the Mach cone 
from the comer. Instead, the procedure must be carried out 
using source lines of quadratic strength. The source strength 
f(x) can then be approximated smoothly so that false pressure 
jumps do not occur. A single source liae of this type repre- 
sents the flow past a slender pointed body with a cusped nose 
(see fig. 12), as is clear from the slender-body approxi- 
mation (equation (90)). 
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METHOD OF (SOLUTION FOE SMOOTH BODIES 

The second-order solution will be described first for bodies 
having continuous slope. Modifications for treating sharp 
cornel’s will be discussed in the next section. 

The procedure is indicated in figure 13. The axis is 
divided into intervals by choosing points with abscissas 
at each of which a source line is to begin. Good accuracy 
is usually obtained if the interval length is not greater than 
/S times the local radius. The tangency condition will be 
imposed on the surface of the body at the points P n , which 
lie on the Mach lines from the points £«. 




Figuhe 12.— Body formed by source line of quadiatlc strength. 


/r-xjp 





For pointed bodies, the first-order solution is started with 
a conical source from the origin which gives the proper 
conical tip. This potential and the derivatives which are 
required are 


<p 0 =— C 0 x (sech -1 1— Vl— f 2 ) <po : 
<Po =“-6' 0 sech“ l f 


Cn 1 


,/i — t- 


<Po„ 


P«r r'- 


* Vl-f 2 
<3 Co 1 
' x t-Jl—t* 

- 1 - 


Y (91) 


where 


x 


v l — 


V 1 — j9 2 « a + €* sech 'jSe 


Here Co is the same as the A of equations (69) to (73), e 
being the tangent of the semivertex angle of the conical 
tip. No such term is required for the semi-infinite body. 

The subsequent procedure is the same for either body. 
Quadratic source lines are started from each of the points 
£i, &, and so forth. For the pointed body & also lies at the 
origin, while for the semi-infinite body it is at — f}a. For 


the nth such source line, the potential and its derivatives 
are given by 


<P*=— C„{x— $«) 2 £(l+h T » S ) sech-'r*— |v/l— 
<pn x = — 2 C a (x - £») (sech - 1 r» — Vl - T, 2 ) 

<Pn=$C n {x- £,) ra T » secll -1 T„^ 

<Pn„= — 2G b scch -1 f» 



where 


The constants C a are determined successively by imposing 
the first- order tangency condition in turn at each of the 
points P n+1 . From equation (80b), the condition is 

at P B+ i (93) 

where the summation begins with n=0 for the pointed body 
and n=l for the semi-infinite body. In this way, values 
of the complete first-order potential and its first and second 
derivatives are calculated at each of the points on the body. 
The velocities due to the particular second-order integral 
can then be calculated at the same points. Differentiating 
equation (85) gives 

\p x —AP [(„ +Nrip r )ip zx -\-‘Px{^i J rNrip xr ) — ‘^ 

► 

'Pt —M * 1 (p+ Nr<p,)<p„ +¥>i[(iV+ 1 )<p, +Mv rr ]~ | <Pr 1 (<Pr +3 r<p„) | 

(94) 

Finally, the second-order complementary function x is 
determined by repeating the procedure used for ip, finding 
new constants such that the second-order tangency condition 
is satisfied. According to equation (19b), the condition is 

lf'r+2xR r =-ff , (l+^i+Sx* I ) &t F»+i (95) 


The second derivatives of x Reed not be calculated. 

The complete second-order perturbation velocities are 
found as the sums of the contributions from \p and x- Then 
the pressure coefficient can be calculated at each point P n 
from equation (26). 

TEEATMENT OF BODIES WITH COBNERS 

Suppose that the meridian curve of the body has a sharp 
comer, which for convenience may be assumed to lie on the 
Mach cone from the origin, as indicated in figure 14. Then 
the method of solution must be modified for two reasons. 
In the first place, the intervals between source lines would 
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have to be chosen extremely small in order to obtain an 
accurate first-order solution behind the comer. This diffi- 
culty can be overcome by adding a new solution which causes 
a sharp deflection of the streamlines. In the second place, 
even if the first-order solution is determined exactly, the 
second-order solution does not yield the Busemann result 
just behind the corner, as it should since the flow is locally 
plane. This defect is remedied by properly canceling a 


r 



FVrt'RE t4.- Body with comer. 


discontinuity which arises in the second-order solution at 
the corner. 

These two modifications require special solutions of the 
first-order equation which along the Mach cone from the 
origin have discontinuities in velocity in the first case, and a 
discontinuity in potential in the second case. Such solutions 
can be found by approximating to equation (83) in the vicin- 
ity of the Mach cone. Imposing the condition that there x 
is only slightly less than fl/f and keeping only leading terms 
in the difference r—&R leads to an Abel integral equation 
for the source strength. Inverting the integral equation 
shows that a potential having discontinuous nth derivatives 
results from a source distribution along the axis which is 
proportional to Setting /(x)=x“ - * in equation (81) 

gives 


£«-5 fift 

&x,r)= — . - = 

Jo A (a— £)2-0’r 2 

i x-pr) a f l _ 


(96) 


'Phis integral represents the analytical continuation of the 
hypergeometric function, so that, except for a constant 
factor, 



Here t is the conical variable introduced in equation (14), 
and K and E are the complete elliptic integrals of the first 
and second kinds with modulus t=[(l— f)/(l+01 v * 

From the tangeney condition (equation (80b)), it can be 
shown that in order to account for the comer the solution 
given by equation (98) should be multiplied by 

(g/-fi /)[ l+( y> )J .... 

jB+i?/ 


Here RJ and /?/ are the slopes of the meridian curve just 
ahead of and behind the comer and (<Px) a is the value of <p c 
ahead of the comer. The first-order solution can thereafter 
be continued as described in the previous section. 

The second difficulty noted before was that the second- 
order solution is found to be incorrect just behind the comer. 
The proper method of treating this difficulty is to solve the 
case when the comer has been slightly rounded and then 
pass to the limit of a sharp comer. However, the following 
simpler procedure is found to give exactly the same result. 

The particular solution ^ calculated from equation (85) 
is discontinuous along the Mach wave springing from the 
comer. If the discontinuity vanished at the comer, the 
solution could subsequently be revised as in the case of 
plane flow (see fig. 5). However, there is a finite jump in ^ 
directly at the comer, which cannot be allowed. Conse- 
quently, the correction potential x must involve an equal 
and opposite jump. A potential having such a discontinuity 
is obtained by setting n = 0 in equation (97). Then 


2 fa fW 

X r'Vl’Vl+t 


K 




Y »+ 1 i- £ 2^r) (97) 

where a is the radius at the comer. The potential is under- 
stood to vanish except within the downstream Mach cone 
from the origin. The hypergeometric functions occurring 
here can all be expressed in terms of complete elliptic 
integrals with real moduli. 

In the first-order solution, a sharp deflection of the stream- 
lines at a comer is produced by adding a multiple of the 
potential which has discontinuous first derivatives. This is 
found by setting n=l, which gives 


Xx= 


= _1 I [ij_ [W'K W\ 

rz\r 1 —t \ 1 +t ^ 


/£J_ (1 F tA 

70 irx\r E k )i 


(100) 


Adding a suitable multiple of this potential cancels the dis- 
continuity in The second-order solution can then be 
continued as described in the preceding section. It can be 
verified that the pressure jump at the comer has then the 
correct second-order value. 

It is instructive to analyze the behavior of a comer from 
another viewpoint. It was pointed out before that the right- 
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hand side of the iteration equation can be considered to 
represent the effects of a known distribution of sources 
throughout the flow field. In the case of a slightly rounded 
corner, this source distribution will be weak except between 
the Mach lines from the comer. As the comer shrinks to a 
point, the source intensity will increase in that region in 
such a way that the total strength remains constant. In 
the limit, the source distribution will behave like a Dirac 
delta function along the Mach line from the comer. The 
particular integral for plane flow (equation (35)) takes 
account of this impulsive function so that the correct solu- 
tion is automatically obtained. In the case of axially sym- 
metric flow, however, it is clear that the particular integral 
given by equation (85) misses the contribution of the im- 
pulse. It is therefore necessary to correct this shortcoming 
by adding the step-function potential given by equation 
( 100 ). 

COMPARISON WITH NUMERICAL SOLUTIONS 

The accuracy of the second-order solution for bodies of 
revolution can be evaluated by comparison with examples 
calculated using the numerical method of characteristics. 

The first body to be considered is a circular-arc ogive of 
1214-caliber radius of curvature followed by a cylinder, 
which has a half angle of 16.26° at the tip. The sccond- 
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Fionas 16.— Pressure distribution on body of revolution with comer at 2.07*. 


order solution has been calculated for this body at a Mach 
number of 3.24. This represents a severe test of the method 
because the Mach angle is then only 10 percent greater 
than the tip cone angle. Internals were chosen such that 
the points P K lay at 0.1, 0.25, 0.5, 1, 2, and 3.5 calibers. 
The pressure distributions calculated by first- and second- 
order theory are compared in figure 15 with the results of 
various computations by the numerical method of character- 
istics. Of the latter, the result obtained from the interpola- 
tion chart given by Rossow in reference 29 is believed to be 
more accurate than the earlier German computations which 
were taken from the summary report of reference 30. Ex- 
cept near the tip, the second-order solution agrees very 
closely with the numerical results. 

The second body to be considered consists of a cone of 10° 
semi-ve"rtex angle followed by a cylinder. The character- 
istics solution for this body at a Mach number of 2.075 has 
been given by Liepmann and Lapin in reference 31. The 
first- and second-order solutions were calculated beyond the 
corner using the modifications discussed in the preceding 
section. Figure 16 shows the shape of the body, the loca- 
tion of source lines, and the pressure distributions calculated 
by first-order theory, second-order theory, and the method 
of characteristics. Again, the second-order results agree well 
with the characteristics solution. 

SERIES EXPANSION WITH RESPECT TO THICKNESS 

An alternative method of solving the exact perturbation 
equation (equation (4)) by successive approximations is to 
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assume that the solution can be expanded in powers of the 
thickness parameter e. Thus the exact perturbation poten- 
tial is written as 

(1) (2) <3> (4) 

$=$e+$« 2 -f$£ , +#« 4 + . . . (101) 


Substituting into equation (4) and equating like powers of e 
yields a sequence of equations 


a) a) 0 ) 

■^+$«-0 2 $==O 


(2) (2) (2) 
<£ ro +$«-/9 2 $*r=2 M 



Cl) (l) (1) (l) CD CD "1 

4^ 'j 4’— . -f" J 


;io2) 


which can be solved in succession. The first is again the 
usual linearized equation. This method was applied to plane 
subsonic flow in references 6 and 12. 

Schmieden and Kawalki first pointed out (reference 14) 
that the power series assumed here does not always exist, 
even for plane flow. In general, terms of the form e’ n ln a e 
appear, beginning, with, ^Ine in the third-order solution for 
plane flow and in the second-order solution for axially sym- 
metric flow. Furthermore, for a body of revolution only 
even powers of e arise. Consequently, it is necessary to 
. assume a more general series of the form 

CD C2) (3) (4) (5) 

e+4>e*+ . . . (103a) 

for plane flow and of the form 

Cl) , C2) , (3) C4) CD , C0) „ 

$= ir-l- e + 't’e 0 fn € + ^« 6 + • • • 

(103b) 

for axially symmetric flow. 

On the basis of this assumption, Broderick has developed 
a second-order solution for supersonic flow past slender 
pointed bodies of revolution (reference 32). The analysis 
is rather lengthy, since the simplification resulting from the 
discovery of a particular integral does not appear. The 
results are limited to shapes for which the cross-sectional 
area is given by an analytic function, or at any rate has its 
first four derivatives small of order e, and the first two con- 
tinuous. This is a severe limitation since, for example, 
the two bodies discussed in the previous section are not 
admissible. 

Broderick’s method yields the slender-body counterpart 
of the present second-order theory. Just as the usual first- 
order slender-body results can be derived as asymptotic 
forms of linearized solutions, so Broderick’s second-order 
slender-body results can be obtained by expanding the 
present second-order solution in powers of t and In t for 
small t and retaining secondary as well as leading terms. 
The logarithmic terms arise from the series 

sech -1 f=frij— -yf 2 — • • (0<«1) (104) 


The expansion will now be carried out for the case of flow 
past a cone. 

It is clear from equation (70) that the constant A in the 
first-order solution (equation (69)) is given approximately by 

( 105 ) 
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Substituting this value into equations (73a) and (73b), 
expanding in powers of t and In t, and imposing the tangency 
condition (equation (71b)) shows that 

S=«*+* 4 [(2Ji*-l)Zn ^-A/W+l)-4]+ • • • (106) 

Then according to equation (73), the velocity perturbations 
on the surface of the cone are 


v—**h-*[*<!*$ 




V 37 2 . 

U~e * ln p e + 


(107a) 

(107b) 


Replacing A r by its value from equation (11), the approxi- 
mate pressure relation of equation (27) gives for the pressure 
coefficient on the surface of the cone 


c,-.’(a»! r i)+ 

0 [ < ‘( h D’] (108) 

This is Broderick's result (reference 32, equation (81)). 

This series is compared in figure 17 with the original form 
of the second-order solution which uses the exact pressure 
relation. For the most slender cone, the expansion in series 
causes only a moderate loss in accuracy. For more practical 
thicknesses, however, the expansion reduces the accuracy to 
such an extent that for the cone of 20° semivertex angle 
Broderick’s solution is inferior to the first-order result (with 
the exact pressure relation). The reason must be that the 
iteration process itself converges more rapidly than do the 
subsequent expansions which axe required to reduce it to the 
slender-body series form. Hence terminating all expansions 
at terms of the order of those retained in the iteration process 
results in an unnecessary loss of accuracy. 

THREE-DIMENSIONAL FLOW 

PABTIAL PARTICULAR INTEGRAL 

It might be hoped that a particular integral, which so 
greatly simplifies the iteration for plane and axially sym- 
metric flows, could be found for the general three-dimensional 
case. The various methods of existing first-order theory 
could then be applied immediately to the problems of second- 
order flow past such shapes as inclined bodies of revolution 
and three-dimensional wings. 

A part of such a particular integral is found at once, being 
common to the two special cases. Consider the three- 
dimensional iteration equation (equation (12)), which may be 
written 


<j> -f- — fP4>xz = AP {2N& l <p x <p xx —20 L <p x <Pzx +2 (pyVzy -f-2^^4- 

( 109 ) 


A=P+ . . . 
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(a) Semi- vertex angle 5°. 
(o) Semi- vertex angle 15°. 


(b) Semi-vertex angle 10°. 
(d) Semf-vertex angle 20 °. 


Figgee 17— Effect of expanding in series upon second-order pressure on a cone. 


It can be readily verified that except for the term in N and 
the triple products (the last three terms) a particular integral 
is given by 

Tp* = MVPl (HO) 

which appears in both equations (35) and (85). 

The iteration equation is thereby reduced to 

4>y> + 4>ti — ^4>xi —M 1 { 2 fPNip z <p xz + PxVw + 2 <p v tp,<p tt + p*V«) 

( 111 ) 


It has not been possible to find a particular integral of lliis 
equation in terms of the first-order potential. The solutions 
for plane and axially symmetric flow do not appear to suggest 
a generalization. On the other hand, there is no assurance 
that such an integral cannot be found. When the triple 
products are negligible, the right-hand side of equation (111) 
vanishes for 7= — 1 (iV=0). However, investigation of the 
previous solutions indicates that the idea of hero taking 
7=— 1 is not legitimate. 
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In the absence of a complete particular integral, the re- 
duced iteration equation (equation (111)) must be attacked 
by more conventional methods. In principle, it is always 
possible to find a particular integral of a linear nonhomoge- 
neous equation with the aid of the fundamental solution 
associated with the differential operator. For the three- 
dimensional wave operator which occurs here, the funda- 
mental solution is 


V(z - £) 2 - j9 2 [(y - tj ) j +(2 - ffl 


which can be interpreted as the potential at any point 
(x, y, z) lying inside the downstream Mach cone from a unit 
supersonic source at (£, ij, f). With the aid of Green’s 
formula, it can be shown that a particular integral of 


is given by 


•IVf ■)- $« — -F(x,y,z) 


(113) 


m x y 2 \ = l 1 ( r 


f) W f 

r,y+( z-ffl 


(114) 


where the integration extends throughout that portion of 
the forward Mach cone from the point (x, y, z) within which 
F is defined. 

In practice, the integration indicated in equation (114) is 
generally not feasible. For example, even the simplification 
of axial symmetry reduces equation (114) only to a double 
integral of F(x,r) multiplied by a complete elliptic integral 
of complicated argument. Avoiding such integrals by dis- 
covery of the particular solution clearly represents a great 
simplification in this case. 8 

In the following sections, one example of a three-dimen- 
sional solution will be given, and the possibility of treating 
other shapes will be discussed thereafter. 


Flow PAST AN INCLINED CONE 

The problem of a cone at an angle of attack illustrates the 
use of separation of variables to reduce the three-dimensional 
iteration equation to tractable form. 

Two alternative coordinate systems are suitable for bodies 
of revolution at an angle of attack. In wind axes the body 
is inclined, while in body axes the stream impinges jsn the 
body obliquely. The latter system is simpler for first-order 
problems and is probably better for the second-order solution 
also. However, wind axes will be used here, since otherwise 
the iteration equations must be rederived. 

To simplify the solution, it will be assumed that the angle 
of attack a is so small that its square can be neglected. This 
will give a solution nonlinear in the body thickness but linear 
in a, and will therefore yield the initial slope of the lift curve 
correct to second order. The coordinate system is indicated 
in figure IS. 6 

* Compart i] r the two methods would lead to the evaluation of definite Integrals Involving 
complete elliptic Integrals, which might be of some Interest. 

* In previous versions of this work (references 1 and J) an oblique transformation was ap- 
plied which effectively unyaws the ails of the body. Howeve, this additional transforma- 
tion has been found to complicate rather than simplify the analysis, and has been omitted 
here. 


To this approximation the surface of the cone is given by 

r=I?(x, 0) = [€—a(l + 6*) cos 0]x (115a) 

f=T(0)=j3[6-a{l + £ i )cosff] (115b) 

The first-order equation for the conical perturbation poten- 
tial is given by the left-hand side of equation (17): 

(l-f*)?«+f+^T=0 (116) 

The solution required here is the sum of potentials for a coni- 
cal line source and dipole (reference 3, p. 74) and has the 
form 

p=— .d^seeh"^— Vl— ^ — — t sech -1 tja cos 0 

(117) 


r 



From equations (16) and (19a) the tangency condition is 
found to be 


=[«-«(! + **) cos 0](l-f5-fp,) &tt=T (118) 

Substituting from equation (117) and expressing values on 
the cone in terms of their values at f=|8e by means of Taylor 
expansions, it is found that 


Vl — j8 2 e*+e s sech. 


C=jS 6 2 (i+0 


(1 +2t I )yi-0 , e*+/SVsech- 1 


(119a) 


(119b) 


According to equation (110), a partial particular integral 
of the second-order iteration equation is, in conical form, 

= il/^-<?,)= Axlf 2 { A sech" 1 <(sech _I f-VT^f 2 )- 

^r^-yl^F sech _i f -1-t* ,.,..1 ) 

(7 3- 2 — ^ f(sech 1 1) 2 acosfl) — 

cos 9 (120) 


There remains to solve the reduced iteration equation given 
by equation (111), which becomes 


(1 -t^„+^+^AM^A(2N ^=^-,S 2 ^ O^)- 


( 121 ) 
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Figure 19.— Comparison of various solutions for normal force on a cone. 


This is reduced to two total differential equations by setting 

?(<, 0)=<£o(f)+$i(t)a= cos 8 (122) 

Here is associated with the axial component of the free 
stream, and ?i with the cross flow. The first of these is 
known from the previous solution for the symmetrical cone 
(equation (72)). The equation for the cross-flow term ?i is 

(1 -tw+if-h IACM< (n ( ^|==+ 

(123) 

Setting 

?!(f)=to(f) (124) 

reduces this to a linear first-order equation in u which can bo 
integrated to find that 

?!=Z> sech -1 tJ+ACM s |^3A T 2 ~ + 

Nt (sech" 1 *)*+§ (125) 

The tangency condition separates into the two conditions 
^W+? 0 , )=«[l+(^o-f^)+(5o-f5o')] at t=pe (126a) 


0 + (^+?i)=^(l+6 J )W'+<Ao")- 

[l+^+^-^o'+WI (12Cb) 

The first of these is the previous relation (equation (71b)) 
which determined the constant B in equation (72). Simi- 
larly, the second of these determines the constant I) in 
equation (125). 

« 

SERIES EXPANSIONS FOR PRESSURE AND NORMAL FORCE ON IN CUN ED 

CONE 

Numerical results have been calculated only for the case 
in which the solution is expanded in powers of t and In t (the 
slender-body approximation). Carrying out the expansion, 
the constant D is found to bo 

D= 2/3« s £l +2M 2 eHn ~~(3APN+^j e ! + . . .J (127) 

Then calculating the velocity components from equation 
(117) and the pressure from equation (27) gives, on the 
surface of the cone, 

<7, =«?,)„- 4e [l-e*(il/ 2 fa~-|il/ 2 -H)]acos0+ . . . 

( 128 ) 
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Here (C f ) 0 is the value for zero angle of attack, given by 
equation (108). Integrating gives the normal-force coeffi- 
cient, based on cross-sectional area: 

forec =2 r ( ipln j.3 J/l+1 \ + "] 

| P .L-*(area) L V 2 > J 

(129) 

This result has been obtained also by Lighthill (reference 33), 
who has calculated the lift on bodies of revolution having 
analytic meridian curves by assuming a series expansion for 
the velocity potential. 

Stone (reference 34; has developed a solution for inclined 
cones which is linearized with respect to a, but otherwise 
exact. Kopat (reference 35) has published tables of the 
numerical results of Stone’s theory. A comparison of equa- 
tion (1291 with this exact theory and with Tsien’s first-order 
solution (reference 36) is shown in figure 19 for 5° and 10° 
cones. The earlier discussion of series expansions suggests 
that the agreement might improve if the solution were not 
expanded in series. 

SHOCK-WAVE POSITION FOB INCLINED CONE I 

Just behind the Mach cone the velocity components are 

-(j?)' 1 = 2:i/2AV ( 1 - 8 ^ aCosfl ) ( 13 °) 

For simplicity, using equations (119), A and (J have here 
been approximated by « a and 20e a . Comparing equations (74) 
and (77 J, it is seen that if these were the velocities just 
behind the shock wave, then the difference between the 
shock-wave angle and the Mach angle would be 

X-sin- 1 i=Md/ 2 A r V(l-8/3a cos 0) (131) 

Hence the ratio of the angular rotation of the shock wave 
to that of the cone would be 

(132a) 

a 

Recently Lightliill has derived a simple expression for the 
shock-wave position for any conical body lying inside 
the Mach cone (reference 25). From his results it is found 
that the preceding relation is incorrect, the correct expression 
being again one and one-half times as large, so that 

-=3( 7 + 1)* e i =l2M 2 8 2 A ri t i (132b) 

a p 

For moderate cone angles this latter result agrees well with 
the exact values calculated by Ivopal (reference 35) from 
Stone’s theory'. Again the discrepancy indicates failure of the 
second-order solution near the Mach cone. 

EXTENSIONS OF THE THEORY 

Two important classes of problems have only' been touched 
upon here. One is wings; the other, inclined bodies of 
revolution. The problem of inclined bodies has recently' 


been studied further in reference 37 . The iteration equations 
were there rederived in body' rather than wind axes, which 
simplifies the tangency condition. Again only a partial 
particular integral could be found. 

The possibility of discovering particular integrals of the 
iteration equation might be investigated more systematically. 
If none can be found for general three-dimensional flow, 
special cases such as conical flow should be studied. 

TREATMENT OF WINGS 

One of the most useful applications of first-order theory' is 
to thin flat wings. For conical wings, the reduced iteration 
problem can be transformed, by' the standard conical theory 
(reference 18), into the problem of solving Poisson’s equation 
inside a circle. This case has recently been studied by 
Moore (reference 38), who has calculated results for the 
nonlifting wing of triangular plan form lying inside the Mach 
cone and of symmetrical wedge section. Although the two 
extremes of the plane airfoil and the slender cone show that 
compressive pressures are reinforced in the second approxi- 
mation ? Moore calculates a reduction for an intermediate 
sweptback case. 

Two difficulties can be anticipated. First, if the wing has 
subsonic edges, infinite velocities arise there, so that the 
assumption of small perturbations is violated. It is known 
that in first-order theory this is no essential objection, 
since the pressure is found correctly' except in the immediate 
neighborhood of the singularity, and the integrated values of 
lift and moment are correct to first order. Selunieden and 
Kawalki (reference 14) and Kaplan (reference 13) have 
indicated that this result extends to the second approximation 
for subsonic flow, so that probably' no real difficulty' exists. 

Secondly, if the wing has supersonic edges, the failure of 
the iteration process along Mach lines from the apex can be 
expected to affect the surface pressures. Again it is possible 
that integrated values will be correct to second order. 
Otherwise, it may' be possible to adjust the solution in those 
regions, in a manner similar to that shown in figure 5. 

HIGHER APPROXIMATIONS 

It seems unlikely that third or higher approximation 
would ever be justified. Other neglected factors, chiefly 
viscosity and heat conduction, should ordinarily be considered 
first. However, the Busemann second-order result has been 
extended to third and even fourth order (reference 39), 
and various writers have considered the third approximation 
for plane subsonic flow (references 7, 9, and 11). If a third 
approximation should be considered worthwhile, the iteration 
could be repeated. Again the eases of flow past a curved 
wall ancl a cone would serve as helpful examples. 

APPLICATION TO SUBSONIC FLOW 

The iteration equation and the particular integrals are in 
no way restricted to supersonic flow. The particular solution 
for plane flow might profitably' be compared with the sub- 
sonic solutions of references 5 to 4. 

Recently Kaplan has shown (reference 40) that the par- 
ticular integral for plane flow (equation (35)) can be derived 
formally', using complex variable theory, and has similarly 
obtained the third-order particular integral. The latter is 


'These relations were given Incorrectly In references l and 2 . 
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not applicable to supersonic flow except in special cases which 
are free of shock waves. 

The particular integral for plane flow has recently been 
applied to several subsonic problems by Harder and Klunker 
(references 41 and 42). 

The particular solution for axially symmetric flow makes 
possible a second-order solution for bodies of revolution at 


subsonic speed. In this case, the integral equation can be 
treated by the methods used for the airship problem. 

Ames Aeronautical Laboratory 
National Advisory Committee for Aeronautics 
Moffett Field, California, April 2, 1952 
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SYMBOLS 


constant reference radius for body of revolu- 
tion 

abscissa at which source distribution for body 
of revolution begins 

constants determined by boundary conditions 
local speed of sound 
constant coefficients of series 
normal-force coefficient 
pressure coefficient 

complete elliptic integral of the second kind 
with modulus &=[(1 — f)/(l+f)] 1/s 
source-distribution functions for body of 
revolution 

known right-hand side of (ra+1) st-order iter- 
ation equation 

continuous function of order unity which 
vanishes for z=0 

arbitrary functions of one variable 
complete elliptic integral of the first kind with 
modulus £=[(1— f)/(l+f)] 1/s 


free-stream Mach number 



7+1 M 2 
2 0 2 


local static pressure 

points on body of revolution at which lan- 
gency condition is imposed 
local speed of flow' 
radius in cylindrical coordinates 
radius of meridian curve of body of revolution 
function defining surface of body 

conical variable 

value of t on surface of body 
perturbation velocity components in Cartes- 
ian or cylindrical coordinates 
free-stream velocity 

Cartesian coordinates with x in free-stream 
direction 

value of y on surface of planar body 
angle of attack 
VM 2 — 1 

adiabatic exponent 

angular rotation of shock wave on cone due 
to angle of attack 

parameter small compared w r itli unity 
azimuthal variable in cylindrical coordinates 
angle of shock w'ave from free-stream direction 



p 

T» 


<t> 

$<*> 

(»> 

X 

'l'* 

a 

Ct> 


(ft) 


0 

1 

CO 

a 

b 

c 


abscissa of origin of nth source line 
local density 

conical variable referred to a=£* as origin 
rather than r=0 

first-order (linearized) perturbation potential, 
same as <f> a) 

second-order perturbation potential, same 
as 

exact pertubation potential 
nth-order perturbation potential 
nth term in series expansion of perturbation 
potential 

complementary function for second-order iter- 
ation equation 

particular integral of second-order iteration 
equation 

partial particular solution for three-dimen- 
sional flow 
(See equation (110).) 
complete velocity potential 
auxiliary variable 
(See equation (124).) 

SUPERSCRIPTS 

result of nth-order solution 
value in conical form 
(See equation (15).) 

SUBSCRIPTS 

associated with axial flow' 
associated with cross flow 
free-stream conditions 
value just ahead of corner 
value just behind comer 
differentiation in direction of cross-wind com- 
ponent of normal to body surface 
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